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Abstract 



We examine an optimization problem in the tropical mathemat- 
ics setting, which is motivated by practical problems in just-in-time 
scheduling. The problem is to minimize a nonlinear objective function 
defined through a multiplicative conjugate transposition on vectors in 
a finite-dimensional semimodule over a general idempotent semifield. 
To solve the problem, we first derive a lower bound for the objective 
function, and then find vectors that make the objective function equal 
to the bound. Under fairly general conditions, an explicit solution is 
given in a closed vector form. The result is applied to solve scheduling 
problems formulated as minimizing a function in the form of the span 
scminorm, subject to linear constraints. As an illustration, numerical 
examples of optimal scheduling are presented. 
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nonlinear objective function, span seminorm, project scheduling. 

MSC (2010): 65K10, 15A80, 65K05, 90C48, 90B35 

1 Introduction 

Tropical (idempotent) mathematics is concerned with the theory and appli- 
cations of semirings with idempotent addition. Tropical mathematics has 
its origin in seminal works [TJ [21 [SJ HI [5] , where it is introduced as a natural 
tool to represent and solve real-world problems in operations research, in- 
cluding scheduling problems examined in [2JE]. Over the past few decades, 
there has been significant progress in the field, which is reflected in several 
monographs and a wide range of research papers, among them most recent 
monographs El EJ Efl 03] . 

Since an early researches [12\ [T3] , optimization problems that are denned 
and solved within the framework of tropical mathematics constitute an im- 
portant research domain in the field. The problems are usually formulated 
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as to minimize (maximize) functions on finite-dimensional semimodules over 
idempotent semifields, subject to constraints in the form of linear equalities 
and inequalities. Both linear and nonlinear objective functions are consid- 
ered. 

One of the objective functions, which is encountered in the problems, 
is the span (range) seminorm, which is given, in the ordinary setting, by 
the maximum deviation between components of a vector. This function 
is used as an optimality criterion for some problems in a range of areas 
from the analysis of Markov decision processes [141 [Td] to the form-error 
measurement in precision metrology |16|. 117] . In the context of tropical 
mathematics, the span seminorm is introduced by |18|ll9j. where it is called 
the range seminorm. 

The span seminorm appears in |2CHI21| in a tropical optimization problem 
coming from machine scheduling. The problem deals with machines, each 
producing a component for final products, and aims at finding starting time 
of production for each machine so that the completion time of all machines 
is spread over a shortest possible period of time. A solution to the problem 
is given in a form that involves both semifields K m ax,+ and K m in,+ • 

In this paper, we consider a more general problem of just-in-time schedul- 
ing |22[ |23| . It is formulated in the common setting of project scheduling 
as having a set of activities operating under various precedence relations 
between their initiation and completion times. The problem is to design a 
schedule that provides, as far as possible, a single common completion time 
of all activities, and thus can be solved by minimizing a span seminorm. 

We represent the precedence relations in terms of the semifield (R max ,+ 
by linear vector equalities and inequalities. The span seminorm is written 
in a straightforward, if not linear, vector form. As a result, we arrive at a 
constrained optimization problem with a nonlinear objective function, which 
involves a conjugate transposition operator, and with linear constraints. 

The above problem is taking as a motivating example and a starting 
point to define and solve a new tropical optimization problem in a rather 
formal setting. We exploit the fact that the application of the solution 
for linear inequalities in [24, 10J allows to reduce the problem to an un- 
constrained problem with new variables. We examine an extended version 
of the unconstrained problem formulated in terms of a general idempotent 
semifield. To solve the latter problem, the solution approach developed in 
[25\ \TU[ [26] is implemented based on the derivation of a sharp lower bound 
for the objective function. 

We give an explicit solution to the extended problem under fairly general 
conditions, which is represented in a compact vector form in terms of the 
carrier semifield. The above scheduling problems are then solved as par- 
ticular cases. Specifically, an alternative direct solution is suggested to the 
machine scheduling problem in [201 121~] . 

The paper is organized as follows. It begins with a motivating problem 
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drawn from just-in-time scheduling in Section [2j Furthermore, we give a 
brief introduction to basic definitions, notation, and preliminary results in 
tropical mathematics in Section [3] to provide a formal framework for subse- 
quent results. Section 0] suggests the main results and includes definition and 
solution of general optimization problems with nonlinear objective functions. 
Application of the results to optimal scheduling problems are presented in 
Section [SJ accompanied with illustrative numerical examples. 

2 Motivating example 

We start with a real- world problem taken from project scheduling and in- 
tended to both motivate and illustrate further results. The problem arises 
in just-in-time manufacturing and aims to design a schedule that minimizes 
the maximum deviation between completion time of activities in a project, 
subject to various activity precedence constraints. For more details and ref- 
erences on project scheduling, and specifically on just-in-time scheduling, 
see, e.g., |22H23]- 

Consider a project which consists of n activities (jobs, tasks) operating 
under start-finish and start-start precedence constraints. The start-finish 
constraints require that a minimal time lag be held between initiation of one 
activity and completion of another. Each activity is assumed to complete as 
early as possible to meet the constraints. The start-start constraints specify 
a minimal time lag between initiation of any two activities. The problem is 
to find a schedule that provides, as far as possible under the constraints, a 
single common completion time for all activities. 

For each activity i = 1, . . . , n, let Xi be an initiation time, yi be a com- 
pletion time, and Cjj be a minimum possible time lag between initiation of 
activity j = 1, . . . ,n and completion of i. Given time lags Cj,- , the comple- 
tion time of activity i must satisfy the start-finish precedence relations 

•Ej Cij — Vii J ' = 1) • • • ; ^> 

with at least one inequality holding as equality. Note that we assume ca > 
for all i. Provided that Cy is not given for some j , we put = — oo. 
Now we combine the relations into one equality of the form 

max (xj + Cij) = yi. 

l<j<n 

Furthermore, let dij be given minimum possible time lag between initi- 
ation of activity j and initiation of i. Once again, we assume dij = — oo if 
no lag is specified for i and j . Due to the start-start constraints, we have 
relations 

•Ej ~\~ dij ^ Xj, j — 1, . . . , 7Z, 
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and then rewrite them as one inequality 



max (xj + dij ) < x,- t . 

l<j<n 

Finally, we define an objective function to formulate the design of op- 
timal schedule as an optimization problem. An optimality criterion, which 
naturally shows if a schedule provides a single common completion time for 
all activities, is the maximum deviation between the completion times, or 
the span seminorm 

max yi — min yi = max yi + max (—yi). 

l<i<n l<i<n l<?<n l<i<n 

Now we can formulate an optimization problem of interest as follows. 
Given a« and bij for all i,j = l,...,n, find x%, . . . ,x n that 

minimize ( max yi + max (—yi)), 

l<i<n l<i<n 

subject to max (xj + cy) = yi, m 

l<j<n ' ' 

max (xj + dij) < Xi, i = 1, . . . , n. 

l<j<n 

Note that a particular case of the problem is examined and solved in 
\20\ [21] in the framework of tropical mathematics. A solution is given in a 
combined form that requires using both semifields K m ax,+ and K m in,+ • 

Below we show how the problems can be represented in a different trop- 
ical mathematics setting and then directly solved in a compact vector form. 



3 Preliminary results 

In this section we give a brief overview of main algebraic definitions, notation 
and preliminary results, which provide a basis for subsequent solution of 
tropical optimization problems and their applications to project scheduling. 

Both concise introductions and thorough presentations of the theory and 
methods of tropical mathematics are given in various forms in a range of 
works, including [H El O El [271 HI]- Below we mainly adhere to results 
in |24[ [TO] , which offer a useful framework for getting direct solutions in a 
compact form. For additional details, one can consult other publications 
listed above. 

3.1 Idempotent semifield 

Let X be a set that is closed under two associative and commutative oper- 
ations, addition © and multiplication ®, and equipped with their neutral 
elements, zero and identity 1. Addition is idempotent, which implies 
that x © x = x for all x € X . Multiplication is distributive over addition 
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and invertible, that is, each x € X + , where X + = X \ {0}, has inverse x _1 
to satisfy x~ l ® x = 1. Since X + forms a group under multiplication, the 
structure (X, 0, 1, ©, <g>) is commonly referred to as the idempotent semifield. 

The integer power is introduced as usual. For any x £ X + and integer 
p > 0, there are defined x° = 1, P = 0, x p = x p_1 <g>:c, and :r~ p = (z _1 ) p . 

In what follows, the multiplication sign ® is dropped for simplicity. The 
power notation is used in the sense of the above definition. 

The idempotent addition produces a partial order, by which x < y if 
and only if x © y = y . The partial order is assumed to extend to a consistent 
total order over X . From here on, the relation symbols and the minimization 
problems are thought in the context of this order. 

As examples of the general semifield under consideration, one can take 
Rmax,+ = (R U{-oo}, -oo,0, max, +), K min ,+ = (KU{+oo},+oo,0,min,+), 
Kmax,x = (K+U{0},0, l,max, x), and IR m i n , x = (K+U{+oo}, +oo, 1, min, x), 
where R is the set of reals, R + = {x € R\x > 0}. 

Specifically, the semifield K m ax,+ has the null = — oo and identity 
1 = 0. Each x 6 R gets its inverse x~ x given by — x in standard notation. 
For any x,y £ R, the power x y is equal to the arithmetic product xy . The 
order, which is induced by addition, corresponds to the natural linear order 
on R. 

3.2 Matrix algebra 

Now we take into consideration matrices over X and denote the set of matri- 
ces with m rows and n columns by X mxn . A matrix with all entries equal 
to is called a zero matrix. A matrix is row (column) regular, if it has 
no zero rows (columns). A matrix is regular, if it is both row and column 
regular. 

For any conforming matrices A = (a,ij) , B = (bij) , C = (cij) , and scalar 
x matrix addition, matrix and scalar multiplication are routinely defined as 

{A © = ciij © bij, {B © C}ij = bik c kji {xA}ij = xciij. 

k 

For any matrix A, its transpose is denoted by A T . 

For any nonzero matrix A = (ay) € X mxn , we introduce a multiplicative 
conjugate transpose A € x raxm with elements a~- = aj* if ciji ^ 0, and 
a~- = otherwise, for all i = 1, . . . , n and j = 1, . . . , m. 

Let us turn to square matrices in X nxn . A matrix that has all diagonal 
entries equal to 1 and off-diagonal entries to is an identity matrix denoted 
by i". For any matrix A, its trace is given by 

n 

tr A = Ojj . 

i=l 
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The matrices that have only one column (row) are routinely referred to 
as column (row) vectors. We denote the set of column vectors of order n by 
X n . 

A vector having all components equal to is a zero vector. A vector is 
regular if it has no zero components. 

Let x be a regular column vector and A be a matrix. It is not difficult 
to see that the vector Ax is regular only when the matrix A is row regular. 
Similarly, the row vector x T A is regular only when A is column regular. 

As usual, a vector y is linearly dependent on vectors xi,..., x rn if there 
are scalars ci,...,c m S X such that y = c\X\ © • • • © c m x m . Specifically, a 
vector y is collinear with x if it holds that y = cx for some scalar c. 

For any nonzero vector x = (xi) £ X™, its multiplicative conjugate 
transpose is a row vector x = (x^), where xj = x^ 1 if x-i 7^ 0, and 
x~^ = otherwise. The next properties of conjugate transposition are easy 
to verify. 

For any regular vectors x and y of the same size, the componentwise 
inequality x < y implies x~ > y~ and vice versa. 

For any nonzero column vector x, it holds that x~x = t. Moreover, if 
the vector x is regular, then xx~ > I . 

3.3 Solution to linear inequality 

Given a matrix A £ X nxm , consider a problem of finding regular solutions 
x € X" that satisfy an inequality 

Ax < x. (2) 

Below we present solutions given to the inequality in [24\ [TO] and written 
here in a more compact equivalent form. 

For each matrix A € X nxn , we introduce a function 

Tr(A) = tr A © • • ■ ©tr A n . 

Furthermore, provided that Tr(A) < 1, a star operator is defined to 
send A to a matrix A* that consists of linearly independent columns of the 
matrix 

J© A© ••• © A" -1 . 

Lemma 1. Let x be the general regular solution of inequality Then 
the following statements hold: 

1. If Tr(A) < 1, then x = A*u for all regular vectors u. 

2. If Tr(A) > 1, then there is no regular solution. 
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4 Optimization problem 

In this section we present our main result that solves an extended problem 
formulated in terms of a general idempotent semifield. We follow the solu- 
tion approach based on deriving sharp bounds on the objective function in 
tropical optimization problems, which is introduced in [25J and then applied 
in a range of studies, including [TOj, [26] . 

Given matrices A, B € X mxn and vectors p,q G X" 1 , the problem is to 
find regular vectors x € X™ that 

minimize q Ax{Bx)~ p. (3) 

The next statement offers a direct solution to the problem. 

Theorem 2. Suppose A is column regular and B is row regular matrices, 
p is nonzero and q is regular vectors. Denote A = {B{q~ A)~)~ p . 

Then the minimum in problem ([3]) is equal to A > and attained at 
any vector 

x = a(q~A)~, a > 0. 

Proof. To verify the statement, we first show that A is a lower bound for 
the objective function in ([3|), and then present vectors x that provide the 
bound. 

With inequality xx~ > I , we write 

q Axx > q A. 

Since the vector q is regular and the matrix A is column regular, both 
sides in the last inequality are also regular and q~ Ax > for any regular 
x. Therefore, we have 

(q~ Ax)~ l x = {q~ Axx~)~ ^ 

Multiplication by B from the left gives 

(q~ AxY 1 Bx < B(q~A)~. 

Considering that the matrix B is row regular, both sides of the inequality 
are regular vectors, and thus 

q Ax(Bx) > (B(q~A)~)~. 

After right multiplication of both sides by the vector p, we finally have 
the lower bound in the form 

q Ax{Bx) p > (B(q~ A)~)~p = A > 0. 

It remains to verify that x = a(q~A)~ yields the bound for any a > 0. 
Indeed, substitution into the objective function and identity x~x = 1 give 

q Ax(Bx) p = q~ A{q~ A)~ (B{q~ A)~)-p = (B(q~A)~)~p = A. □ 
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We conclude this section with solutions of two particular cases of the 
above problem. First assume that A = B = I and p = q = 1, where 1 
denotes a vector having all components equal to 1 . We arrive at a problem 
that is to 

minimize 1 xx 1. 

Application of Theorem [2] immediately gives A = 1 as the minimum in 
the problem, which is attained at any vector x = at for all a > 0. 

Finally, we examine a problem that underlies the design of optimal sched- 
ules to be given below. We put B = A, p = q = 1 , and consider a problem 

minimize t T Ax{Ax)-t. (4) 

With Theorem [21 we immediately obtain the following result. 

Corollary 3. Suppose A is a regular matrix and denote A = (A(t T A)~)~t . 

Then the minimum in problem is equal to A and attained at any 
vector 

x = a(l T A)~, a > 0. 

5 Applications to optimal scheduling 

We are now in a position to put the scheduling problem described above 
into the framework of tropical mathematics and then give a direct solution 
to the problem in a compact vector form. 

5.1 Representation of scheduling problem 

Consider problem (fTJ and note that its representation in ordinary notation 
involves only operations max, addition, and additive inversion. On this 
basis, it is natural to represent the problem in terms of the semifield K m ax,+ • 
First we represent constraints in terms of K m ax,+ as scalar equalities and 
inequalities 

n n 

C-i j %j — Ui j dij Xj ^Xj, i — 1 , . . . , 71. 

3=1 3=1 

With matrix-vector notation 

C = (cij), D = (dij), x = (xi), y = {yi), 
the scalar constraints take the form 

Cx = y, Dx < x. 
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Furthermore, we rewrite the objective function in (pQ). Since, for (R max ,+ , 
it holds that 1 = (0, . . . , 0) T , the objective function is readily given by 

Finally, by combining the objective function with the constraints, we 
arrive at a problem 

minimize 1 yy 1, 

subject to Cx = y, (5) 
Dx < x. 

5.2 Solution of scheduling problem 

A direct solution to problem ([5]) under general conditions is given as follows. 

Theorem 4. Suppose that C is a regular matrix and D is a matrix with 
Tr(.D) < 1. Denote A = (CD*(l T CD*)-)-±. 

Then the minimum in problem ([5]) is equal to A and attained at 

x = aD*(l T CD*y, y = aCD*(l T CD*)-, a > 0. 

Proof. It follows from Lemma Q] that the inequality constraints in problem 
© have the solution x = D*u for all vectors it > of appropriate size. 

With the solution, the equality constraints become y = CD*u with 
u > 0. 

Substitution of y in the objective function of ([5]) leads to problem ([3]) 
with A = B = CD* , p = q = 1, and an unknown regular vector u. 

Since D* > I, the matrix CD* is regular. Application of Corol- 
lary gives the minimum A = (CD*(1 T CD*)~)~1 attained at u = 
a(l T CD*)-. 

It remains to go back to the original vectors and get x = aD* (1 T C D*)~ 
and y = aCD*(l T CD*)~ for all a > 0. □ 



Finally consider a problem from \20\ [21] , that can now be represented as 

(6) 



minimize 1 yy 1, 
subject to Cx = y. 



As a direct consequence of the above results, we get the following solu- 
tion. 

Corollary 5. Suppose C is a regular matrix and denote A = (C(1 T C)~)~1 . 
Then the minimum in problem ^ is equal to A and attained at 

x = a(l T C)~, y = aC(l T Cy, a ><D. 
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5.3 Numerical examples 

To illustrate the results obtained, we consider an example project of three 
activities under constraints given by matrices 

/ 0-2 1 
D = 2 

\ -l oo 

where the symbol = — oo is used for ease of exposition. 

First we do not take into account the start-start constraints so as to 
solve reduced problem (|6j). After calculating the vectors 

/ -4 \ / 

1 T C = ( 4 3 3 ) , (1 T C)- = I -3 j , C(1 T C)- = I 

we apply Corollary [5] and immediately arrive at a solution 

(~*\ (°\ 

A = 0, x = a -3 , y = a , a>0. 





Note that in this situation we really get a just-in-time schedule with a 
single common completion time of all activities. 

Let us now incorporate the start-start constraints given by D into the 
problem. We take D to calculate 

\ / -1 -2 

D 2 = [ 1 -2 1 I , D 3 = [ 0-1 2 ) , Tr(.D) = 0, 




and then evaluate the sum 



I D D 2 = | 1 




-1 

Since the first and the third columns in the obtained matrix are collinear, 
we drop the last column to define the matrix D* . We successively get 




D* = | 1 | , CD* = | 4 3 | , 1 J CD* 

Furthermore, we have vectors 

D*(1 T CD*)- = [ -3 ) , CD*{t T CD*)- = [ 
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Application of Theorem 2] gives the results 




a > 0. 



The solution offers an schedule that is optimal with respect to the span 
seminorm. Note, however, that the constraints in the problem give no way 
for the schedule to provide a single common completion time of all activities. 
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